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' Abstract 

o : 

04 ' A model is introduced to describe guided propagation of a linear or nonlinear pulse which encounters a localized 
5h ' nonlinear defect, that may be either static or breather-like one. The model with the static defect directly applies to 
an optical pulse in a long fiber link with an inserted additional section of a nonlinear fiber. A local breather which 
gives rise to the nonlinear defect affecting the propagation of a narrow optical pulse is possible in a molecular chain. 
J — i In the case when the host waveguide is linear, the pulse has a Gaussian shape. In that case, an immediate result of 
■ its interaction with the nonlinear defect can be found in an exact analytical form, amounting to transformation of the 
incoming Gaussian into an infinite array of overlapping Gaussian pulses. Further evolution of the array in the linear 
host medium is found numerically by means of the Fourier transform. An important ingredient of the linear medium 
is the third-order dispersion, that eventually splits the array into individual pulses. If the host medium is nonlinear, 
the input pulse is naturally taken as a fundamental soliton. The soliton is found to be much more resistant to the 
Ph action of the nonlinear defect than the Gaussian pulse in the linear host medium, for either relative sign of the bulk 
and local nonlinearities. In this case, the third-order-dispersion splits the soliton proper and wavepackets generated 
by the action of the defect. 

> 

in 

O 

o 

O 

> 

X 



< 



1 



I. INTRODUCTION 



The interaction of traveling solitary pulses (which, in particular but not necessarily, may be solitons, that we realize 
here as pulses in nonlinear media maintaining a steady shape while propagating) with various local imperfections or 
pinned dynamical excitations is a problem of fundamental importance. For direct experimental observation, the most 
straightforward case is the propagation of a pulse in an optical fiber (linear or nonlinear) , in which a strong localized 
nonlinear defect represents an inserted piece of a dispersion- shifted fiber (DSF) [Q, that has finite nonlinearity and 
negligible dispersion. This configuration, which may be quite useful for optical telecommunications [Q, is described 
by the following version of the nonlinear Schrodinger (NLS) equation, 

iUz - -(32UrT + -^PsUttt + 7 1^1^ u + TS (z) \uf M = 0, (1) 

2 o 

where m(z,t) is a local amplitude of the electromagnetic wave, z is the propagation coordinate, 

T = t- z/Vo (2) 

is the so-called reduced time, Vq being the group velocity of the carrier wave, f32 and P3 are coefficients of the second- 
order and third-order group-velocity dispersion (GVD) in the fiber, and 7 is the nonlinearity coefficient of the host 
(system) fiber. The strength of the localized nonlinear defect is F = 7dsf-^^, where 7dsf and L are the nonlinearity 
coefficient and actual length of the above-mentioned finite- length DSF inserted into the host fiber at the point z = 0; 
this insertion may be represented by the delta- function in Eq. (|l|), as the DSF's dispersion is negligible Note 
that the nonlinear coefficient 7 is always positive, as it is induced by the Kerr effect in the optical fiber, that always 
has the sign corresponding to self-focusing ||]. On the other hand, the GVD coefficient /32 may be either positive or 
negative, which corresponds, respectively, to normal and anomalous dispersion The third-order dispersion (TOD) 
coefficient (i^ is frequently neglected, unless (32 is very small or if the pulse is very narrow in the r-domain. However, 
TOD will play an essential role in the present work. Note that is usually positive in optical fibers 

A qualitatively different version of the model containing the local defect pertains to the case when it describes 
the interaction of the propagating modes (pulses) with a localized defect in the form of a pinned breather . These 
are spatially localized, time-periodic dynamical states, which are ubiquitous in nonlinear physical systems ranging 
from quasi-one dimensional polymers and charge-density- wave materials (e.g., metal- halogen electronic chains 
[0) to Josephson ladders Scattering of a moving solitary pulse on a breather is an important process in many 
physical, chemical and biological systems which combine charge, spin and energy localization with transport of these 
quantities by pulses. Detailed understanding of this type of the scattering will not only yield valuable information 
on the dynamical properties of breathers, such as their mobility and stability, but also reveal the dynamical response 
of materials in which the breathers are excited. In particular, the interaction of a narrow optical pulse with a pre- 
existing electron-phonon breather in a lattice chain is a process of considerable fundamental and technological interest, 
since it may enhance optical nonlinearities of the material, and thus the efficiency of the second- and third-harmonic 
generation. Indeed, there is some evidence of such an effect in conjugated polymers such as polyenes 

Another relevant example of the breathers is provided by the rotational dynamics of certain chemical groups in a 
molecular chain ]To| . Analysis of the pulse scattering on them, in conjunction with recent developments in ultrafast 
(femtosecond) spectroscopy and chemistry , will enable an efficient use of the rich photophysics of functional 
optical and electronic materials. Damage tracks in certain mica minerals and sputtering on crystal surfaces have 
also been attributed to moving breathers In addition to mixed- valence transition- metal complexes 0, there 

is experimental evidence for localized breather- like states in antiferromagnetic chains p^ , as well as in the above- 
mentioned Josephson junction arrays and ladders ||^. 

The one-dimensional (ID) model introduced and considered in this work is a first step in a systematic study of the 
pulse-breather scattering, the eventual objective being to elucidate the interaction of pulses and breathers in 2D and 
3D nonlinear dynamical lattices. However, the formulation and consideration of the corresponding multidimensional 
models is a complicated issue, which is beyond the scope of this work. 

The simplest generalization of Eq. (|l|) for the case when the local nonlinearity is induced by a small-size breather 
oscillating at a frequency lo is 

iuz P2UTT + —PsUttt + 7 1""!^ u + rS (z) cos^ (uit) ■ \u\^ u = 0. (3) 

2 6 

Equation (^) implies that oscillations of the breather modulate the strength of the corresponding local nonlinear 
defect, but do not change its sign, which is a natural assumption in the case of the cubic nonlinearity. Note that, 
despite the difference between the time t and the reduced time r, see Eq. (^, it makes no difference what definition 
of time is used in the argument of cos^ in Eq. (||), as r = i at z = 0. 
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Both equations and (0) conserve the norm of the wave field (which has the physical meaning of energy in the 
application to fiber optics 0), 

/ + 00 
|u(r)|^dT. (4) 
-oo 

Additionally, the model (|l|) including the static defect conserves the field momentum, 

/+00 
UTU*dT, (5) 

the asterisk standing for the complex conjugation. 

The simplest and perhaps most fundamental version of both models and is the one with all the nonlinearity 
concentrated at the point z = 0, while the bulk nonlinearity is negligible, i.e., 7 = 0. It appears, rather surprisingly, 
that, unlike the more technically involved case of the fully nonlinear model with 7 > 0, which was studied in a part in 
Ref . [0 , in the case 7 = the interaction of the moving pulse with the nonlinear defect has not been yet considered. 
Therefore, the detailed study of the pulse-defect interaction in the model with 7 = is the first objective of the 
present work. Then, we will also consider the full nonlinear model with 7^0, concluding that, as a matter of fact, 
the most interesting results can be obtained just in the model with the linear host medium, 7 = 0. 

In the linear medium, it is natural to take an incident pulse in the form of a Gaussian, which is an eigenmode of 
the linear system. The problem of the interaction of a Gaussian with the delta-like nonlinear defect can be partially 
solved in an analytical form, which is done below in section 2. In fact, the analytical result can be obtained in an 
advanced form for the case of the static defect with uj — 0, which corresponds to Eq. (|^) and, as it was explained 
above, has direct application to fiber optics. A basic result produced by the analytical consideration is that the initial 
Gaussian pulse, passing through the nonlinear defect, generates an infinite series of strongly overlapped Gaussian 
pulses. At this stage of the analysis, TOD becomes a crucially important ingredient of the model: we show that it 
lends each Gaussian pulses its own velocity, which will bring about eventual separation of the pulses. In section 2, 
the velocity generated by TOD is found in an analytical form, which yields an exact result for an initial stage of the 
evolution. The splitting of the solution into an array of pulses is considered in detail in section 3. 

In section 4 we display results of the fully numerical solution for the case when the incident Gaussian pulse interacts 
with the dynamical defect [the breather, see Eq. (^)]. In this case (w ^ 0), the eventual pattern is much less regular, 
consisting of separated pulses with random shapes, rather than Gaussian-like ones. However, if the frequency uj is 
very large, one should expect that the dynamical defect described by the last term in Eq. (||) may be replaced by its 
averaged static counterpart. In accordance with this expectation, the numerical computations demonstrate return to 
a regular pattern for very large values of uj. 

In section 5 we present a solution for the most general model, which includes the bulk nonlinearity at z ^ 0, i.e., 
7 7^ in Eq. ( ^). In this case, it is natural to take the incident pulse as a soliton of the corresponding homogeneous 
NLS equation, rather than a Gaussian pulse, and the interaction may only be simulated numerically. Results of 
simulations of the interaction of the soliton with the defect turn out to be quite different from those in the case when 
the host medium was linear: the soliton is found to be much more resistant to the action of the local nonlinear defect 
than the Gaussian pulse in the linear medium. Unless the defect is extremely strong, the pulse remains virtually 
intact, the TOD term separating it from small-amplitude wave packets generated by the local defect. However, in 
the absence of the TOD term, the effect of the perturbation may accumulate and destroy the soliton. Thus, the 
third-order dispersion plays a crucially important role in both versions of the present model, i.e., for the linear and 
nonlinear host medium. 

If the defect is very strong, it may split the soliton into two. The results obtained for the model with the nonlinear 
host medium are not sensitive to the value of the defect's intrinsic frequency w, being essentially the same for the 
static defect and its dynamical counterpart. Moreover, the soliton is found to be stable against the action of the local 
defect irrespective of the relative sign of the bulk and local nonlinearities. These results, obtained for the case when 
the input signal is a soliton, rather than a linear pulse, is another manifestation of the well-known general principle, 
according to which solitons are very robust eigenmodes of nonlinear media. 



II. ANALYTICAL CONSIDERATION OF THE INTERACTION BETWEEN A GAUSSIAN PULSE AND 

A NONLINEAR DEFECT. 
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A. The pulse 



We start the consideration with the simplest version of the model, viz., Eqs. (|l|) or (|) with 7 = /Jg = 0. At 
z 0, we thus have a linear Schrodinger equation, which gives rise to the well-known exact solution in the form of a 
Gaussian pulse (which is termed a coherent state in quantum mechanics), 



u{z,t) 



exp 



c(l + 2iP2Cz) 2 
l + {2p2czf 



(6) 



where c > and Aq are arbitrary real constants. The constant c determines the width of the Gaussian pulse and its 
chirp (i.e., the imaginary part of the coefficient in front of in the argument of the exponential, which shows a slope 
of the local frequency across the pulse |^). 

Note that the solution (||) has zero velocity in the present reference frame [which is defined by Eq. (||)]. An exact 
solution for a moving pulse can be generated from Eq. (^ by the action of a boost (Galilean transform) 



u(z, t) u(z,T — sz) exp I —T-z - 



IS 



P2 



(7) 



where s is a real velocity-shift parameter. It is important to notice that, for the boosted pulse, the values of its 
velocity, momentum (||) and energy (Q) are related, irrespective of the particular form of the pulse, in a simple way. 



-^pulse 



pulse 



(8) 



B. Pulse acceleration by the third-order dispersion and numerical verification 

Before proceeding to the consideration of the passage of the pulse (||) through the nonlinear defect, we need to 
understand how a free pulse will move under the action of the TOD term in Eqs. (^ or ( The consideration of this 
issue will help to understand how an array of pulses splits in the presence of the TOD term. To this end, a solution is 
sought for by means of the Fourier transform, which yields the following integral representation for it: 



u{z,t) 



Ao 



exp 



iujT + i 

4c 
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du 



(9) 



[setting z = 0, this expression goes over into Eq. (^ taken a,t z — 0]. In the limit of large values of z, the integral (||) 
is dominated by a contribution from a vicinity of the stationary-phase point, which is 



Wo 



T 



ZT 



J2 ^ 

In particular, a contribution of the stationary-phase point ( p^ to the phase of the pulse is 



(^,r) 



2/32Z 6/3|z2- 



(10) 



(11) 



Direct simulations show that, under the action of TOD, the initial Gaussian pulse is gradually destroyed, generating 
a long "tail" , while a part of the wave packet may still be interpreted as a surviving pulse. Using the expression (|lT 
one can find a relation between the momentum [see Eq. (E)] and energy of the pulse- like part of the wave packet. 



P - 

-'pulse — ^ -L^pulse- 
^P2 



(12) 



Note that, as the net momentum of the wave field must be conserved, a "recoil" momentum — Ppuiso is carried away 
by the above-mentioned tail. The comparison of Eqs. (^ and ( p^ ) yields an analytical prediction for the asymptotic 
value of the velocity shift acquired by the surviving pulse under the action of the TOD term in the limit z — > 00: 



s = (l/2)/33C. 



(13) 
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To this end, a natural definition of the position of the wave-packet's center of mass is adopted, 



T, = E-^ T\u{T)\'dT, (14) 

^ — OO 

where E is the net energy defined by Eq. (|^), and the derivative drddz may be regarded as a velocity of the pulse. 
In the case of the free propagation, one can derive an exact relation 

where P and E are the net momentum and energy of the wave field. Comparing this to Eq. (||), we conclude that, 
in the absence of TOD, the velocity ( |l5|) is identical to the boost parameter (in particular, dTc/dz = if P = 0); 
however, this identity is broken by TOD, that is why acceleration of the pulse by TOD is observed. If the TOD term 
is treated as a small perturbation, one can analytically calculate the second term on the right-hand side of Eq. ([T^), 
using the expression (|), which is an exact solution in the absence of TOD. As a result, we find the velocity of the 
pulse in the case when the net field momentum vanishes, P = 0, which is true for the initial Gaussian pulse (^: 

[note that this expression does not contains neither z nor [32, despite the fact the solution (|^), used for the calculation 
of the right-hand side of Eq. (|l^), does depend on z and /32]. The expression (|l6|) shows that the velocity lent to the 
pulse by TOD linearly depends on the pulse's parameter c, hence initially overlapping pulses with different values of 
c [see Eq. ( pl| ) below] are expected to separate under the action of TOD. 

The analytical result (|2^) was checked against direct numerical simulations of the linear version of Eq. (|l|). For 
instance, in the case = 1, /^s = 0.1, and c = 1, it was found that the pulse was indeed moving at a constant velocity, 
the value of which exactly coincided with that given by Eq. (|l^), in the interval < z < 25. At larger values of z, 
the absolute value of the velocity decreases, which can be explained by the fact that the TOD term essentially alters 
the shape of the pulse at that later stage of the evolution. 

C. Passage of the pulse through the nonlinear defect 

Our next aim is to consider transformation of the Gaussian pulse passing the nonlinear defect. Obviously, in an 
infinitesimal vicinity of the defect (as \z\ 0), only the first and last terms should be kept in Eqs. ( |l|) and (|^), which 
yields a simplified equation, 

^ ^ iT5{z) cos' {lot) ■\u\^u. (17) 
az 

To solve Eq. ( [l7|) , we represent the solution as u{z) = a{z) exp [i0(z)], with real amplitude a and phase (f). Substituting 
this into Eq. (^7|), one immediately finds that da/dz — 0, and 

^ =r5(z)cos2 {ujT)-a\ 
oz 

A solution to the latter equation is obvious, 

0(z = +0, r) - (j){z = -0, t) = To' (t) ■ cos' (wr) , (18) 

where we take into regard that a may be a function of r. 

Thus, we take the input pulse at the point z = —0 in the general form [cf. the expression (|^)], 

u{z = -0) = Aa exp [- (co + ibo) t'] , (19) 

where cq > determines the initial width of the pulse, and bo is its initial chirp. The substitution of the expression 
( p^ ) into Eq. (|lj) yields the form of the pulse appearing after the passage of the nonlinear defect: 

u{z = +0, t) = Ao exp [- (co + ibo) t' + iTAl exp (-2cot^) • cos^ (cjt)] . (20) 
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Further analytical consideration for the general case of the dynamical defect (breather), with w 7^ 0, is extremely 
cumbersome. Therefore, in the rest of this section and in the next one, we focus on the static case, to = 0. The aim 
will be to realize a result of the further evolution of the transformed pulse ( pO| ) , governed by the linear equation (|l|) 
Including the TOD term) at z > 0. To this end, we notice that the expression (20) can be expanded into an infinite 
series: 



+00 



u{z = +0) = ^0 exp (-ibor^) ^ (n!)"^ (iFA^)" exp [- (1 + 2n) cqt^] 



(21) 



n=0 



Comparing it to the exact fundamental-pulse solution (|6|), we conclude that the expression (|2^), if considered as an 
initial condition to the linear equation (|^) with 7 = /Ja = 0, gives rise to a superposition of an infinite number of 
Gaussians with the values of the width constant c„ = (1 + 2n) cq, and with the common initial value 60 of the chirp. 
The evolution of the wave packet ( pT| ) can be presented in a relatively simple form in the case when the initial chirp 
is absent, 60 = (and TOD is neglected, /^a = 0): 



+00 



uiz, t) = Ao exp {-iboT^) J2 , ^ exp 

nly'l - 2z/32Co (1 + 2n) z 



Cq {1 + 2iP2Co {I + 2n) z) _ 
1 + (1 + 2nf {2f32CQzf 



(22) 



[in the case w 7^ 0, Eq. (|2^) shows that the Gaussian in each term of the initial series ( |2l| ) with n 7^ is additionally 
multiplied by \cos{ujt)]^" , which will make the subsequent result much more complex than that given by Eq. ( |2^)]. 

Thus, Eq. ( p^ ) gives an exact solution to the nonlinear model equation (|l|) in the case when the incident pulse 
has no chirp and 7 = /^s = 0. However, while taking the input pulse to be chirpless, and disregarding the bulk 
nonlinearity (7 = 0)are quite acceptable assumptions, the TOD term may not be neglected, as, without this term, 
the exact solution ( p2|) remains strongly degenerate. Indeed, centers of all the Gaussian pulses, the superposition of 
which constitutes this solution, exactly coincide, staying at r = [note that this degeneracy is not lifted if the input 
pulse has nonzero chirp, nor if the bulk nonlinearity (7 7^ 0) is added]. On the other hand, Eq. ( |l6|) shows that the 
TOD term will lend each pulse its own velocity, depending on the initial width parameter c of the pulse. Obviously, 
this will eventually split the superposition (^2|) of the overlapping Gaussians into an array of separating pulses. 

However, the analytical prediction (^6|) for the TOD-induced velocity shift of each pulse is only valid for finite 
values of z. Therefore, to find an actual shape of the evolving wave train, it is necessary to solve numerically the 
linear Schrodinger equation with the TOD term. 



1 i 

lUz - -p2Urr + -^PzUttt = 0, 

2 



(23) 



with the initial condition in the form (20). Results generated by Eq. (|23| ) are presented in the next section. 



III. GENERATION OF THE WAVE TRAIN BY THE THIRD-ORDER DISPERSION IN THE CASE OF 

THE STATIC NONLINEAR DEFECT 



As Eq. (|2^) is linear and has constant coefficients, it can be solved by Fourier transform. Thus, the numerical part 
of the solution amounts to the computation of the Fourier transform 



ua{uj) 



+ OC 



exp(icL'r)w(z — +0, r)dr 



for the initial configuration (|2C|), and subsequent computation of the inverse Fourier transform 

/ + 00 
exp(— ia;T)u(z, ijj)dT, 
-oo 

where, as it immediately follows from Eq. (p3|), 

u(z, uj) = ua{uj) exp i (^j32Uj'^z - i/?3W^z 

In Fig. 1 we display the profiles of |u(t)|, computed at the points z = 10, 20, 30, and 40 for a case when the static 
nonlinear defect is weak, having F = 0.1. The figure demonstrates that, in accordance with the analysis presented in 
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the previous section, the initial wave packet tends to spht into an array of regular Gaussian-like pulses. The splitting 
actually takes place for larger values of the nonlinear-defect's strength, as is shown in Fig. 2, which pertains to F = 1. 
For the same case, the eventual shape of the pulse array is shown in more detail in Fig. 3. 

The result for a still stronger nonlinear defect, with F = 10, is displayed in Fig. 4. In this case, the splitting into 
pulses takes a violent character, with appearance of huge gradients and formation of a very sharp front, the latter 
effect being accounted for by the interplay of the strong local nonlinearity and TOD. These results may be explained 
by the fact that, as follows from Eq. (pT|), the number Timax of the largest-amplitude pulse in the series grows ~ F 
with increase of F, hence the width W of an individual pulse decreases ^ l/'s/F, which gives rise to the large gradients 
revealed by the subsequent evolution. We also note that, although the TOD coefficient /Ja remains small, the size of 
the TOD term in Eqs. @) and (|) grows as with the decrease of the pulse's width. This explains the formation 

of the abrupt front under the action of the asymmetric TOD term. 



IV. GENERATION OF THE PULSE ARRAY BY A LOCALIZED BREATHER 



An example of the splitting of the initial Gaussian as a result of its interaction with the dynamical defect (with the 
same strength F = 1 as in the case shown in Fig. 2, and the frequency ui = 1/2) is shown in Fig. 5. The splitting is 
displayed in more detail by blowups collected in Fig. 6. 

Comparison of Figs. 2 and 5 shows that the static defect and the dynamical one with a moderate value of the 
frequency produce quite similar results. Taking larger values of the frequency strongly changes the situation: as is 
shown in Fig. 7 and in the accompanying blowup (Fig. 8), the same value of the defect strength as in the cases 
displayed in Figs. 2 and 5, i.e., F = 1, but combined with uj — 10, gives rise to an essentially more disordered pattern. 
Note that this pattern is disordered in a way essentially different from that observed as a result of the action of a 
strong static defect (lo = 0) , cf. Fig. 4. In particular, strong asymmetry of the pattern and sharp fronts are not 
found in the present case, and the local gradients are not as huge as in Fig. 4. The absence of those features in the 
present case is easy to understand, as they may only be generated by a large value of F, as explained above. 

Keeping to increase a; at a fixed value of F, we have concluded that the situation shown in Figs. 7 and 8 is rather 
similar to that observed at F = 1 and oj — 100 (not shown here). On the other hand, it is obvious that, if the 
frequency is extremely large, one should be able to replace cos^ (ujt) in Eq. (^ by its average value 1/2, thus reverting 
to Eq. (|l|) for the static nonlinear defect. To verify this argument, in Fig. 9 we present the results of the numerical 
computations for F = 1 and oj = 1000, which are supported by the blowup shown in Fig. 10. From these pictures, it 
is evident that, in the case of extremely large lu, the situation is indeed nearly the same as in the case of the static 
nonlinear defect, cf. Figs. 2 and 3. 



V. TRANSFORMATION OF A SOLITON BY THE LOCAL DEFECT IN THE NONLINEAR HOST 

MEDIUM 

In all the cases considered above, the input pulse was taken in the form of a Gaussian, since the host medium in which 
this pulse propagated before the collision with the nonlinear defect was linear, the Gaussian being its eigenmode. The 
situation is completely different in the case when the host waveguide is itself nonlinear, as in that case the incoming 
pulse must be a sech soliton [|l|. The action of the localized nonlinearity on the soliton is described by the same 
general expression ( p^ ) as above; however, the result cannot be interpreted in such a straightforward manner (even 
if w = 0) as it was done above for the case of the input Gaussian pulse, see Eq. (^). In fact, in the case of the 
nonlinear host medium all the analysis following the use of expression (^8|) for the defect-induced phase change of the 
input soliton can only be performed numerically. 

As is commonly known, in the homogeneous part of the nonlinear model (^) or (^ (at z ^ 0), the soliton may only 
exist if the second-order dispersion is anomalous while the nonlinearity is self-focusing, or vice versa i.e., only if 
/327 < . In the simulations, we fixed (32 = +1, and took different negative values of 7, see below. The input soliton 
was always taken as a fundamental one with a fixed width corresponding to these values of the parameters: 



m(z = -0,t) = |7|-i/^sechT. (24) 

If the TOD term is neglected {Ps = 0), the nonlinear defect produces a strong perturbation around the soliton. 
However, the TOD term helps to separate the soliton and the perturbation, as is illustrated by Fig. 11 for the case 
/Jg — 0.3. This figure shows a large r-domain, in order to demonstrate the evolution over a large propagation distance; 
as is seen, the perturbation spreads out indefinitely, while the soliton remains essentially intact. 
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The comparison of Fig. 11 with Fig. 2, which shows the interaction of the Gaussian pulse in the hncar host 
medium with the local nonlinear defect that has the same strength, F = 1, suggests a conclusion that the soliton in 
the nonlinear host medium is much more resistant to the action of the nonlinear defect than the Gaussian pulse in 
the linear host medium. In fact, this conclusion is strongly supported by many other simulations with different values 
of the parameters. 

Another difference of the present case from that for the Gaussian pulse in the linear host medium is that the result 
of the action of the nonlinear defect on the soliton is less sensitive to the defect's intrinsic frequency. For instance, 
if the defect is dynamical with u) = 1, while the other parameters take the same values as in the case shown in 
Fig. 11, the evolution (not shown here) seems nearly the same as in Fig. 11. with a difference that the generated 
perturbation takes a somewhat larger portion of energy from the soliton (as a result, the soliton reappears after the 
collision with the amplitude « 0.8, to be compared with the amplitude « 1 in Fig. 11). Nevertheless, in this case 
too, the perturbation remains rather small and it does not strongly affect the soliton. 

With further increase of the dynamical defect's frequency, the perturbation again becomes smaller, and the general 
picture is reverting to that corresponding to the static defect with lu — 0. This trend can be easily explained by the 
self-averaging of the dynamical defect similar to that observed above in the case of the linear host medium. However, 
unlike that case, where the self-averaging was evident only for extremely large values of the frequency, u ~ 1000 (see 
Fig. 9) , in the present case simulations (not shown here) clearly show that a; = 10 is sufficient for the self-averaging 
to manifest itself. 

The above results were obtained for the cases with F > and 7 < 0, which implies that the localized nonlinearity 
is self-dcfocusing if the bulk nonlinearity is self-focusing, or vice versa. It is also interesting to consider the case when 
both nonlinearities have the same sign, i.e., F < 0. An example is displayed in Fig. 12 for F = —1 (this example 
pertains to the dynamical defect with w = 1, but the results are virtually the same for the static defect with w = 0). 
In this case, the action of the defect produces a weaker effect, as the soliton reappears after the interaction with the 
amplitude « 1, to be compared with the above-mentioned value w 0.8 of the soliton's amplitude found for the same 
values of parameters but F = -|-1. This difference is quite natural, as the local defect which has the same sign of the 
nonlinearity as in the host medium is acting to compress the soliton stronger, while the local defect with the opposite 
sign of F was acting to stretch the soliton, which produces a more destructive effect, helping some radiation (wave 
packets) to escape from the weakened soliton pulse. 

Taking a much stronger nonlinear defect (with F = 11, see Fig. 13), we conclude that, naturally, the soliton loses 
a larger part of its energy to the generation of the perturbation separating from it. Nevertheless, the soliton survives 
even in this case (cf. Fig. 7, which shows a complete and fast chaotization of the wave field in the case of the collision 
of the Gaussian pulse with a strong defect, having F = 10, in the linear host medium). In the cases when F is large 
but negative (not shown here), when the nonlinearity of the defect has the same sign as in the host medium, the 
soliton's losses are smaller than in the case displayed in Fig. 13, which is quite similar to what is shown in Fig. 12 
for a moderately strong defect with F < 0. 

VI. CONCLUSION 

Wc have introduced a model describing collision of a pulse in a linear or nonlinear waveguide with a strong nonlinear 
local defect, that may be either static or breather-like. The model with the static defect should directly apply to an 
optical pulse in a long fiber-optic link with an inserted section of a nonlinear (dispersion-shifted) fiber of an arbitrary 
length. On the other hand, a local breather, which gives rise to the nonlinear defect affecting the propagation of 
narrow optical pulses, may be realized in a molecular chain with electron-phonon coupling. 

In the case when the host waveguide is linear, the pulse was naturally taken as a Gaussian. A result of its 
interaction with the nonlinear defect was found analytically, amounting to its transformation into a "lump" consisting 
of an infinite number of overlapping Gaussian pulses. Further evolution of the lump in the linear medium is generated 
by the corresponding Fourier transform. An important ingredient of the medium is the third-order dispersion, that 
splits the lump into an array of individual pulses; the velocity shift lent to an initially Gaussian pulse by the TOD term 
was found in an analytical form. The influence of the intrinsic frequency, in the case when the defect is a breather, 
was also investigated. 

The full numerical solution for the model in which the host medium is nonlinear, and the input pulse is taken as a 
fundamental soliton shaped by this medium, produces a result quite different from that for the linear host medium: 
the soliton is much more resistant to the action of the local nonlinearity. If the local defect is not very strong, the 
soliton remains essentially intact, the third-order-dispersion separating the soliton and small wave packets generated 
by the collision. Even if the local nonlinearity is very strong, the soliton survives, losing a limited part of its energy. 

Beyond straightforward application of the results reported above to fiber-optic links, the model introduced in this 



8 



work, and current developments in experimental techniques, such as ultrafast spectroscopy [|ll|,y_2| of optical and 
electronic materials , the results presented here may help to understand the dynamics of scattering of pulses on 
static and breather- like defects in many physical systems, e.g., conjugated polymers Josephson ladders [||, 

and coupled electron- vibron lattice systems ||l^ . 
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FIGURE CAPTIONS 



Fig. 1. The evolution of the field \u{t)\ after the interaction of the input Gaussian pulse with the nonlinear defect, 
in the case a; = 7 = 0, /32 = 1, /^s = 0.1, and T ~ 0.1. 

Fig. 2. The same as in Fig. 1 for a stronger nonlinear defect, with F = 1. 

Fig. 3. Blowups of the field pattern from the last panel in Fig. 2, clearly showing the formation of an array of 
separated pulses. 

Fig. 4. The same as in Fig. 3 for a very strong static nonlinear defect, with F = 10. 

Fig. 5. The evolution of the field \u{t)\ after the interaction of the input Gaussian with the nonlinear dynamical 
defect, in the case cu = 1/2 , 7 = 0, /32 = 1, /^s = 0.1, and F = 1. Except for uj, these parameters are the same as in 
the case shown in Fig. 2. 

Fig. 6. Blowups of segments of the panels from Fig. 5, showing the gradual splitting of the wave packet into an 
array of pulses. 

Fig. 7. The same as in Fig. 5, but with uj = 10. 
Fig. 8. The same as in Fig. 6, but with uj — 10. 
Fig. 9. The same as in Fig. 7, but with ui = 1000. 
Fig. 10. The same as in Fig. 6, but with a; = 1000. 

Fig. 11. The evolution of the field |m(t)| over a very large propagation distance in a large temporal domain after 



the interaction of the input soliton (EJ) with the nonlinear static defect, in the case 7 — — /32 = —1, P3 = 0.3, and 



F = 1, cj = 0. It is evident that the local defect initially creates a large perturbation around the soliton; however, 
the perturbation is diffused away under the combined action of the second- and third-order dispersions, and a robust 
soliton reappears. 

Fig. 12. The same as in Fig. 11, but with uj = 1 and F = —1. 

Fig. 13. The same as in Fig. 12, but with /^s = 0.5 and F = +11. 
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